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We investigate a compact source of entanglement. This device is composed of a pair of linearly 
coupled nonlinear waveguides operating by means of degenerate parametric downconversion. For 
the vacuum state at the input the generalized squeeze variance and logarithmic negativity are 
used to quantify the amount of nonclassicality and entanglement of output beams. Squeezing and 
entanglement generation for various dynamical regimes of the device are discussed. 
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I. INTRODUCTION 

Quantum entanglement and its consequences puzzled 
physicists since the famous paper by Einstein, Podolsky 
and Rosen jj] . It was recognized in recent years that the 
entanglement is not only an intriguing and mind bog- 
gling feature of quantum mechanics but also a crucial 
and extremely useful resource for information process- 
ing. The various protocols relying on the entanglement 
include quantum teleportation 0], entanglement swap- 
ping Q , dense coding 13 , quantum cryptography Q and 
quantum computing 6]. The efficiency of quantum in- 
formation processing significantly depends on the degree 
of entanglement of the quantum state shared by the par- 
ties involved in a given protocol. It is therefore highly 
desirable to establish reliable sources of pure strongly 
entangled states. 

Quantum optics provides a natural and convenient 
framework for the generation of entangled states, their 
manipulation and measurement. The above-mentioned 
protocols have been originally established for the quan- 
tum systems with two-dimensional Hilbert spaces - the 
qubits. The polarization entangled pairs of photons gen- 
erated by means of the spontaneous parametric down- 
conversion proved to be a very good source of entangled 
qubits and they were employed in a large number of ex- 
perimental demonstrations of quantum information pro- 
cessing. 

Recently, however, a significant attention has been 
paid to the quantum information processing (QIP) with 
systems whose Hilbert space is infinite dimensional. A 
typical example of such a system is a single mode of the 
optical field whose Hilbert space is spanned by the infi- 
nite (but countable) number of the Fock states \n). It 
is in the spirit of this latter approach that one considers 
the quadrature operators x, p and one speaks about QIP 
with continuous variables (CV). Similarly as in the case of 
qubits, the parametric downconversion provides a source 
of CV entanglement. Here, the relevant entangled state 
is the so-called two-mode squeezed vacuum (sometimes 
also called the twin beam) that can be either prepared 
by nondegenerate downconversion or via mixing of two 
single-mode squeezed states on a balanced beam splitter. 
This last-mentioned approach has been in fact used in the 
experiment on the teleportation of continuous variables 



@ . The whole setup for the generation of the CV entan- 
gled state is certainly nontrivial and quite complicated. 
A natural question arises whether it could be possible 
to integrate the generation of the single mode squeezed 
states and the mixing on a beam splitter into a single 
small and compact device that would serve as a source 
of the entanglement. 

In this paper we show that the nonlinear optical cou- 
plers that have been thoroughly investigated during re- 
cent years both theoretically and experimentally (see 
0, |j| and references therein) can serve exactly for that 
purpose. We theoretically investigate the generation of 
squeezing and entanglement in the nonlinear optical cou- 
pler consisting of a pair of nonlinear waveguides oper- 
ating by means of the degenerate parametric downcon- 
version. The waveguides are placed close together so 
that the two modes can overlap and exchange energy via 
linear coupling. We study the dependence of the gen- 
eralized squeezing as well as CV entanglement on the 
coupler's parameters. The entanglement is quantified by 
the logarithmic negativity [Io| . which is an easily com- 
putable entanglement monotone. Finally, the optimum 
phase matching condition that yields the maximal entan- 
glement is analyzed. 



II. THE MODEL AND THE EQUATIONS OF 
MOTION 



The device under study is schematically depicted in 
Fig. 1. The nonlinear coupler consists of two linearly 
coupled waveguides where the degenerate parametric 
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FIG. 1: Sketch of nonlinear coupler formed from two nonlinear 
waveguides A and B with susceptibilities x' 2 '- 
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downconversion takes place. A 'blue' photon with fre- 
quency 2lo may split into a pair of identical 'red' pho- 
tons and also an inverse process may occur. Assuming a 
strong coherent pumping of the modes with frequency 2uj 
(second-harmonic modes) and a linear coupling between 
the modes with frequency lu (fundamental modes) in the 
lossless and nondispersive media, we can describe such a 
device by the following effective interaction momentum 
operator \H\ 
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= n (g* A A% 



g* B B 2 F + 2g L A F B F 



(1) 



where ip^it) is annihilation (creation) operator of the 

fundamental mode in the first waveguide and Bp (Bp) 
is annihilation (creation) operator of the fundamental 
mode in the second waveguide. The constant gA(gB) is a 
product of the strong coherent amplitude of the second- 
harmonic mode in the first (second) waveguide and the 
nonlinear coupling constant, which is proportional to the 

(2) (2) 

quadratic susceptibility x\ (Xk )■ The constants gA 
and gs determine the efficiency of the downconversion 
processes. Linear coupling constant denoted as gi, de- 
scribes the energy exchange between the waveguides by 
means of coupling via the evanescent waves. The symbol 
h denotes the reduced Planck constant and h.c. repre- 
sents Hermitian conjugate terms. 

We shall work in the interaction picture where the 
spatial evolution of an operator O is governed by 
the Heisenberg-like equation (see 0, @ and references 
therein) 
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Here z is a spatial coordinate along the direction of prop- 
agation, Gj nt is an interaction momentum operator and 
[, ] denotes a commutator. Using the momentum opera- 
tor Q , we obtain the Heisenberg equations of motion in 
the interaction picture of the form: 



dAp 
dz 

dBp 
dz 



= 2ig A A F + 2ig* L B F . 
= 2ig B B F + 2ig L A F 



(3) 



For our purposes it is convenient to deal with the quadra- 
ture components x, p rather than with the creation and 
annihilation operators. The quadratures can be directly 
measured via balanced homodyne detection. Recall that 
the quadrature components can be expressed as linear 
combinations of the creation and annihilation operators 
as follows: 
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(4) 



We further express the (generally complex) coupling con- 
stants in terms of their amplitudes and phases, 



g i = \gj\exp(iipj) for j = L, A, B . 



(•5) 



With these definitions at hand, we can derive from 
Eq. © the Heisenberg equations for the quadrature op- 
erators. We write them down in the compact matrix 
form, 
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sin (ipj ) , 


Cj = 


2\gA cos (ipj) . 



(6) 
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Eq. 10 represents a system of linear differential equa- 
tions with constant coefficients that can be solved very 
simply by determining the eigenvalues and eigenvectors 
of the matrix M. Explicit analytical expressions for the 
solution can be found in 

El 

Here we write the solution 

in a formal matrix notation 



= S(*)£(0) , S(z) = e 



Mz 



(9) 



This evolution represents a linear canonical transfor- 
mation of the quadrature operators. All possible lin- 
ear canonical transformations of N modes form a group 
of symplectic (also called Bogolyubov) transformations 
Sp(2N 1 R). In our case, for each z the matrix S(z) is an 
element from the four dimensional representation of the 
symplectic group Sp(4,R). 

One of the most important features of the symplectic 
transformations is that they transform Gaussian states 
again onto Gaussian states. The Gaussian states are 
those whose Wigner function is a Gaussian. Such states 
are thus fully characterized by the first and second mo- 
ments of the quadrature operators. Moreover, the first 
moments (mean values) of the quadratures can be always 
set to zero via appropriate displacements applied locally 
to each mode. Since we are interested in entanglement 
properties of the state generated in the coupler, we can 
focus solely on the evolution of the second moments. 



III. THE GENERALIZED SQUEEZE VARIANCE 

It is convenient to arrange the second moments into 
the covariance matrix V whose elements are defined as 
follows [12: 



(10) 



3 



Here A£j = £j — and = Tr(p^j). In our analysis 
of the coupler we shall assume that all processes are spon- 
taneous, i.e., both fundamental modes are initially in the 
vacuum state. The covariance matrix corresponding to 
this input state is proportional to the identity matrix, 



V(0) = 



1 



10 
10 



2 I 1 
1 



(11) 



The evolution of the covariance matrix corresponding to 
the symplectic transformation 10 is given by the follow- 
ing formula: 



VO) = S(z)V(0)S T (^ 



(12) 



where T stands for the transposition. 

Quantum entanglement is one of the manifestations 
of the nonclassical features of the quantum states. The 
CV quantum state is said to have a classical analogue 
if its density matrix can be written as a convex mix- 
ture of projectors on coherent states. Such states cannot 
be entangled and they can exhibit only classical correla- 
tions. Since the nonclassicality is a necessary prerequisite 
for the entanglement 01 it is of interest to investigate 
whether and how much nonclassical are the states gen- 
erated in the nonlinear coupler. This issue has been ad- 
dressed in many recent papers devoted to the analysis of 
the quantum-statistical properties of light propagating 
in nonlinear couplers 0, . Various kinds of nonclassical 
behaviour have been predicted, such as the generation of 
light with sub-Poissonian photon number statistics, and 
single- mode or two- mode squeezed light. 

For Gaussian states, the nonclassicality criterion is par- 
ticularly simple because these states are nonclassical if 
and only if (iff) they are squeezed. More formally, we 
say that the (generally TV-mode) Gaussian state is non- 
classical iff there exists a quadrature X that is a linear 
combination of the quadratures Xj and pj of the N in- 
volved modes such that the variance of X is below the 
coherent-state level 1/2. The smallest variance obtained 
as minimum over the variances of all possible linear com- 
binations of the quadratures Xi and pj is called the gen- 
eralized squeeze variance A 12]. This variance can be 
calculated as the lowest eigenvalue of the covariance ma- 
trix o 



X(z) = min{eig[V(z)]} . 



(13) 



The Gaussian state is nonclassical and squeezed iff A < 
1/2. In our discussion of the coupler's behavior we distin- 
guish two main regimes of operation. In the first regime, 
the linear coupling dominates over the nonlinear interac- 
tions 



2|sl| > \9A+g B \ 



(14) 



and the coupler operates below the threshold. If the op- 
posite inequality holds then the nonlinear interactions in 




FIG. 2: The dependence of the generalized squeeze variance 
A on the length z of the symmetric coupler operating below 
the threshold, \gi\ = 2, \sa\ = |<?b| = 0.2. The three curves 
correspond to the phase differences A(j> — (solid line), A<j> = 
tt/2 (dashed line), and A<j> = n (dotted line). 



the waveguides become dominant and the coupler op- 
erates above the threshold. The threshold condition 
(2|<?.l| = \gA + 9b\) was determined by analyzing the 
eigenvalues of the matrix M (see also Ref. [l4|)- These 
eigenvalues are complex below the threshold, which leads 
to oscillatory dynamics typical for the linear coupling. 
Above the threshold, all eigenvalues are real and the dy- 
namics of the coupler is reminiscent of a pure amplifica- 
tion process. 

The behavior of the coupler also significantly depends 
on the phase matching condition. It can be inferred from 
the analytical expressions for the symplectic transforma- 
tion S(z) that there is a single effective phase difference 



A0 = <j>A - <j>B + 20l 



(15) 



that controls the dynamics of the coupler. 

In figure 2 we plot the dependence of the generalized 
squeeze variance A on the length z of the symmetric cou- 
pler (\gA\ = |<7b|) that operates below the threshold. The 
three curves correspond to three different choices of the 
phase difference A(f>. For A<fi = we observe a periodic 
evolution of X(z). In this specific case, the eigenvalues 
of M are purely imaginary, the dynamics is purely os- 
cillatory, and the quantum state of the two modes pe- 
riodically returns to the initial vacuum state. For the 
two other choices A</> = 7r/2 and Aip — tt the parametric 
amplification sets on and the squeeze variance exponen- 
tially decreases with growing length of the coupler. The 
squeezing is fastest for A<f> = ir. For comparison, we plot 
in figure 3 the function A(z) for symmetric coupler op- 
erating above the threshold. It follows that in this case 
the phase difference almost does not influence the dy- 
namics of squeezing and the three curves shown in Fig. 3 
corresponding to three different phase shifts A(f> almost 
coincide. 

For the asymmetric coupler (\§a\ ^ \9b\) below the 
threshold, the eigenvalues of M have nonzero both imag- 



FIG. 3: The dependence of the generalized squeeze variance 
A on the length z of the symmetric coupler operating above 
the threshold, \g L \ = 0.15, \g A \ = \g B \ = 0.2. The three 
curves correspond to the phase differences A(j> = (solid line), 
A<f> = 7r/2 (dashed line), and A(f> — n (dotted line). 



FIG. 4: The dependence of the logarithmic negativity En 
on the length z of the symmetric coupler operating below 
the threshold, \gi\ = 2, \g A \ = |ffs| = 0.2. The three curves 
correspond to the phase differences A(j> = (solid line), A<fi = 
7r/2 (dashed line), and Acj> = n (dotted line). 



inary and real parts giving rise to periodically modulated 
exponential decrease of X(z). The modulations smooth 
out as Acf> goes to ir. Above the threshold X(z) exhibits 
qualitatively the same behavior as in the symmetric case. 



IV. SEPARABILITY AND ENTANGLEMENT 

Having discussed the squeezing of light in the coupler, 
we now turn our attention to the generation of entangled 
states. A state pab of two subsystems A and B is en- 
tangled iff pab cannot be written as a convex mixture of 
product states 

PAB ^ ^TjPjPAj ® PB, 3 , Pj > 0, (16) 
3 

where paj and f>B,j are states of subsystems A and B. 
Since we assume that the fundamental modes are ini- 
tially in pure vacuum states and the evolution is unitary, 
the two- mode state is pure for all z. It holds that every 
pure bipartite state \i/j)ab that is not a product state 
is entangled. While it is easy to decide whether a pure 
state is separable or not, this problem is much more dif- 
ficult for the general mixed bipartite state. Several sep- 
arability criteria have been proposed in the literature. 
Among them, the most powerful and important is the 
Peres-Horodecki (PH) criterion 0, 0, 0] . According 
to this criterion, a quantum state pab is entangled if the 
partially transposed density matrix 

{PAB) T „t a , n p = (PAB)na,mf3 (17) 

is nonpositive, i.e., it has at least one negative eigen- 
value. In general, this is only a sufficient condition on 
entanglement and there exist entangled states with pos- 
itive partially transposed density matrix (the so-called 



bound entangled states). However, in case of two- mode 
bipartite Gaussian states this criterion turns out to be 
both sufficient and necessary 0, 0| . The PH criterion 
is qualitative, it tells us whether the state is entangled or 
separable, but it does not quantify the amount of the en- 
tanglement present in that state. For this purpose, vari- 
ous entanglement measures have been proposed. The PH 
criterion suggests that a function of the negative eigen- 
values of the partially transposed density matrix could be 
a good entanglement measure. Indeed, such measure can 
be constructed and it is called the logarithmic negativity 

m 

E N (p)=\og 2 [l + 2Af(p)]. (18) 

Here the negativity Af(p) is the sum of the absolute values 
of the negative eigenvalues pj of the partially transposed 
matrix Pab 

M{p)=Y J \P3\- (19) 

3 

The logarithmic negativity En has the great advantage 
that it can be easily computed for any bipartite Gaussian 
state. In fact, analytical formula for En has been derived 
for the general mixed two-mode Gaussian state. The 
details can be found in |Toj . here we only mention the 
main results and formulas. 

It is convenient to decompose the covariance matrix as 
follows 

v =(c-b)' ( 2 °) 

where A and B are the covariance matrices of the modes 
on Alice's and Bob's sides, respectively, and C contains 
the intermodal correlations. We must calculate the sym- 
plectic spectrum (c±,C2) of V that can be obtained by 
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FIG. 5: The dependence of the logarithmic negativity En 
on the length z of the symmetric coupler operating above 
the threshold, \g L \ = 0.15, \g A \ = \g B \ = 0.2. The three 
curves correspond to the phase differences = (solid line) , 
A<j> = 7r/2 (dashed line), and Acf> — -k (dotted line). 

solving the biquadratic equation (only positive roots are 
considered) 

C 4 + (dctA + dctB - 2dctC) C 2 + dctV = . (21) 

The logarithmic negativity is a function of the symplectic 
eigenvalues 

2 

E N =Y,F(cj), (22) 
where the function F(c) reads 

m = \ , °, 9 J 0T ^t (23) 

w | -log 2 (2c) for2c<l. y ' 

The evolution of the logarithmic negativity is plotted 
in Fig. 4. The parameters are the same as in Fig. 2, so 
the coupler is below the threshold. As could have been 
expected, for A(j> = we get a purely oscillatory evo- 
lution and En periodically attains its maximum. The 
oscillations are clearly present also for A(j> — ir/2 but 
the general trend is that En grows with z. This is in 
agreement with the fact that for this detuning the state 
becomes more and more squeezed with growing z, cf. 
Fig. 2. Remarkably, En{z) — for all z when Atfi = ir. 
For this particular phase difference, the coupler behaves 
like two decoupled single mode squeezers that do not in- 
teract at all [lj|. Figure 5 shows the behavior of En(z) 
when the coupler operates above the threshold. We can 
see that the oscillatory evolution is replaced with a mono- 
tonic growth of En{z) with z. However, the amount of 
entanglement significantly depends on A(f> even when the 
coupler is above the threshold. This is in marked differ- 
ence with the generalized squeeze variance that is to a 
large extent insensitive to the choice of Acj>, see Fig. 3. 

In the case of the asymmetric coupler above the thresh- 
old En grows monotonously with growing z. For small 
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FIG. 6: The dependence of the optimal phase difference A(j> op t 
and the maximum of the logarithmic negativity En on the 
length 2 of the symmetric coupler operating below the thresh- 
old \g L \=2, \g A \ = \g B \ =0.2. 



z the rapidity of the growth decreases as A<f> goes from 
to 7T. Below the threshold En behaves similarly as in the 
symmetric case. The only difference occurs for A<p = -k 
when En changes periodically and no decoupling of the 
waveguides takes place. 

The results depicted in Figs. 4 and 5 indicate that 
for each z and a given set of parameters, we may tune 
the phase difference such as to maximize the entangle- 
ment. We have performed numerical calculations and 
determined the optimal phase difference A</> opt and the 
corresponding maximal achievable En- The results for 
the coupler operating below the threshold are presented 
in Fig. 6. The A(j) opt exhibits an interesting behavior. 
Initially, it is optimal to set A(f> opt = and this is opti- 
mal up to certain length zq. At this distance, the optimal 
phase difference becomes to deviate from and A(/> opt 
oscillates with z. The figure 6 suggests that for large z, 
A^ p t approaches some fixed asymptotic value. We note 
that similar results were obtained also for the coupler 
above the threshold. 
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V. CONCLUSIONS 

In this paper we have explored the controllable com- 
pact source of squeezing and entanglement formed by 
a nonlinear optical coupler composed of two nonlinear 
waveguides operating by means of spontaneous degener- 
ate parametric downconversion. Firstly, the behavior of 
the generalized squeeze variance in the below the thresh- 
old and above the threshold regimes has been analyzed. 
Secondly, by calculating the logarithmic negativity the 
possibility of the entanglement generation in the coupler 
has been demonstrated. In particular, it was shown that 
the amount of the entanglement can be controlled via the 
effective phase difference <|15|) . Finally, with the help of 



the numerical calculation the optimal effective phase dif- 
ference providing maximum entanglement at the output 
of the coupler has been found. Our theoretical analysis 
thus clearly illustrates the potential utility of the coupler 
for controlled generation of CV entanglement. 
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